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I. INTRODUCTION 



Dynamical chiral symmetry breaking (DCSB) and 
confinement are two fundamental features of Quantum 
Chromodynamics (QCD). It is generally believed that 
with increasing temperature and baryon number density 
strongly interacting matter will undergo a phase transi- 
tion from the hadronic matter to the quark-gluon plasma 
(QGP) which is expected to appear in the ultrarelativis- 
tic heavy ion collisions. These two phases are generally 
referred to as the Nambu-Goldstone phase which is char- 
acterized by DCSB and confinement of dressed quarks 
and the Wigner phase corresponding to QGP in which 
chiral symmetry is partially restored and quarks are not 
confined. Theoretically, these two phases are described 
by two different solutions, the Nambu-Goldstone solu- 
tion and the Wigner solution of the quark gap equation. 
The existence of these two solutions in the chiral limit 
(the current quark mass m — Q) has been shown in the 
framework of Dyson-Schwinger equation (DSE) approach 
of QCD (see, for example, [ll, 0])- However, it is a gen- 
eral conclusion in the previous literature that when the 
current quark mass m is nonzero, the quark gap equation 
has only one solution which corresponds to the Nambu- 
Goldstone phase and the solution corresponding to the 
Wigner phase does not exist [3, |j|- This conclusion is 
hard to understand and one will naturally ask why the 
Wigner solution of the quark gap equation only exists 
in the chiral limit while does not exist at finite current 
quark mass. Furthermore, this conclusion is in fact not 
compatible with the current study of chiral phase tran- 
sition of QCD. In order to see this more clearly, let us 



have a look at the generally accepted QCD phase dia- 
gram (see, for example. Fig. 3 of Ref. [5j). As is shown 
in the QCD phase diagram, it is generally believed that 
at low temperature and high baryon number density the 
chiral phase transition of QCD is of first order and the co- 
existence of the Nambu-Goldstone phase and the Wigner 
phase should appear. It is well-known that in the real 
world the current quark mass is nonzero. If one cannot 
find the Wigner solution of the quark gap equation in 
the case of nonzero current quark mass, this will mean 
that we cannot talk about the co-existence of these two 
phases. This is obviously an unsettled and important 
problem in the study of QCD phase transitions. The au- 
thors of Ref. [6] first discussed this problem and asked 
whether the quark gap equation has a Wigner solution 
in the case of nonzero current quark mass. Subsequently, 
the authors of Refs. [^-Q further investigated the prob- 
lem of possible multi-solutions of the quark gap equation. 
However, as far as we know, this problem has not been 
solved satisfactorily in the literature. In the present pa- 
per we try to propose a new approach to investigate this 
problem. 

The main motivation of the present work is to study 
the Wigner solution of the quark gap equation at finite 
current quark mass and provide a new viewpoint on par- 
tial restoration of chiral symmetry at finite chemical po- 
tential. This paper is organized as follows: in Sect. II 
we analyze the reason why in the previous literature the 
Wigner solution of the quark gap equation does not ex- 
ist in the case of the finite current quark mass and pro- 
pose a new approach to discuss this question. In Sect. 
Ill, based on such an approach, we show in the frame- 
work of Nambu-Jona-Lasinio (NJL) model that the quark 



DSE has a Wigner solution at finite current quark mass. 
Then, in Sect IV we generahze this approach to the case 
of finite chemical potential to study partial restoration of 
chiral symmetry and compare our results with the cor- 
responding ones in previous literature. The results are 
summarized in Sect. V. 



II. QUARK GAP EQUATION AND ITS 
SOLUTIONS 

In order to illustrate our new approach more clearly, let 
us now briefly recall the usual arguments which exclude 
the existence of the Wigner solution of the quark gap 
equation when m ^ Q. The quark DSE under rainbow 
approximation reads as following (in the present paper 
we will always work in Euclidean space, and take the 
number of flavours, Nf = 2 and the number of colours, 
TV, = 3) 



G-\p) = Gi,\p) + 



(2^ 



(1) 

where G{p) is the dressed quark propagator, Go{p) — 
(i^ ■p + m)~^ is the free quark propagator, g is the strong 
coupling constant and Dfj,^{q) is the effective dressed 
gluon propagator. According to Lorentz structure anal- 
ysis, one has 



G-\p)^zM{p^)+B{p'), 

where A{p'^) and B{p^) are scalar functions of p^. 
stituting Eqs. ([2]) into Eq. ([1]), one has 
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where Landau gauge has been employed. From Eqs. ^ 
and (JH one can find when ra — Q there are two dis- 
tinct solutions for B{p'^). One solution is B{p^) ^ 
which describes the Nambu phase, and the other one is 
BijP') = which describes the Wigner phase. However, 
when m ^ 0, it can be easily seen that B{p^) = is 
not a solution of Eqs. ([3]) and ^. From this observa- 
tion one often concludes that when m ^ 0, the quark 
DSE has only one solution corresponding to the Nambu 
phase and the Wigner solution does not exist. Here, it 
should be noted that in obtaining this conclusion one 
has assumed that the dressed gluon propagators in these 
two phases are the same. However, since the features 
of these two phases are so different, it is reasonable to 
expect that the behavior of the dressed gluon propaga- 
tor should be different in these two phases (for example, 
in the familiar liquid-solid phase transition of water, the 



effective interactions between molecules are different in 
the two phases). To see this more clearly, let us look at 
the graphical representation of the DSE for the dressed 
gluon propagator given in Fig. [T] From Fig. [T]it can be 
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FIG. 1. The DSE for the dressed gluon propagator 

seen that the quark propagator can affect the gluon prop- 
agator through quark-loop insertions. Therefore, in prin- 
ciple, since the quark propagators in Nambu phase and 
Wigner phase are quite different, one naturally expects 
that the gluon propagators in these two phases should be 
different, too. Here we would like to stress that this ob- 
servation is model independent. Besides, this observation 
has been verified in the study of quantum electrodynam- 
ics in 2 + 1 dimensions (QED3) by using the coupled DSE 
for the fermion and photon propagators with a range of 
fermion-photon vertices |10| (QED3 has many features 
similar to QCD, such as spontaneous chiral symmetry 
breaking in the massless fermion limit and confinement. 
Due to these reasons it can serve as a toy model of QCD). 
This indicates that one should choose different forms of 
gluon propagator as input to solve the quark propagators 
in the two different phases. Now, the key problem is how 
to choose appropriate model gluon propagators as input 
to calculate the dressed quark propagator in the Nambu 
phase and the Wigner phase, respectively. 

From Fig. [T]it can be seen that one can formally split 
the full gluon propagator into two parts as following 

D,M^Dl^{q) + D'^,M), (5) 

where 1?^^ is the pure Yang-Mills part which includes 
all diagrams without quark loop insertions (which is usu- 
ally called quenched gluon propagator in lattice QCD) 
and D^^, is the quark-affected part which includes all di- 
agrams with quark loop insertions. Obviously, the pure 
Yang-Mills part in Wigner phase should be same as that 
in Nambu phase, whereas in principle the quark-affected 
parts in these two phases should be different. At present 
it is impossible to calculate the two parts D^^{q) and 
D^^{q) from first principle of QCD. So one has to resort 
to various nonperturbative QCD models to express them 
phenomenologically. 



Over the past few years, considerable progress has been 
made in the framework of the QCD sum rule [il|, which 
provides a successful description of various nonperturba- 
tive aspects of strong interaction physics. We naturally 
expect that it might provide some useful clue to the study 
of the model gluon propagator. From the QCD sum rule 
approach the lowest-order contribution of quark conden- 
sate to the gluon propagator is [12] 
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where {ip{y)ip{z)) is the non-local quark condensate and 
(ipip) is the ordinary two-quark condensate, the ellipsis 

2 

represents terms of higher orders in ^ which we ne- 
glect in the present work since we limit our discussion to 
two light flavors u and d. It is evident that the value of 
quark condensate {iptp) is different in the Nambu phase 
and Wigner phase. This makes the gluon propagators in 
these two phases be different. Therefore, in the following 
calculation we can phenomenologically identify A^j,(p) 
in Eq. © as a good approximation of the D'^^{q) part 
in Eq. ^. 



III. NJL-LIKE MODEL AND TWO DISTINCT 
SOLUTIONS AT ZERO CHEMICAL POTENTIAL 

Now, we should specify a model framework to calculate 
the quark propagators in Nambu phase and the Wigner 
phase. The dressed quark propagators are the most el- 
ementary of the n-point Green functions of QCD. It is 
evident that the Dyson-Schwinger equations (DSEs) are 
the natural tool for investigating it in the continuum. In 
particular, it has been shown that DSEs are capable to 
describe the chiral phase transition and deconfinement 
phase transition at finite temperatures and chemical po- 
tential [2, ll3l - [l7J . However, as is shown in Ref. [18l |. 
the Nambu- Jona-Lasinio (NJL) model can capture the 
main physical features of QCD at finite temperature and 
chemical potential. For example, partial restoration of 
chiral symmetry, the critical end point and color super- 
conductivity are all first studied in the framework of the 
NJL model. This is the reason why the NJL model is 
the most widely used QCD model in the study of QCD 
phase transition at finite temperature and chemical po- 
tential (although this model has two defects, namely, it 
can neither accommodate confinement nor is renormal- 
izable). Therefore, as a first step, for simplicity in this 
paper we shall employ the NJL model to study the quark 
propagators in Nambu phase and the Wigner phase. 

In the normal NJL model the following model gluon 



^'D^u{p-q) = Sfj 



■m 



-e{K^ 



(7) 



is employed to calculate the quark propagator, where Mq 
is some effective gluon mass scale and A serves as a cutoff 
and is set to be 1.015 GeV in Ref. [H. This model 
gluon propagator concentrates on the infrared region of 
the interaction which is believed to be vital for DCSB 
of QCD. With such a model gluon propagator Eq. ^ 
becomes 
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The solution of Eq. © is A{p'^) = 1 and B{p'^) = M with 
M being a constant satisfying the following equation 
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i?i(M',A2), 



(9) 



where Di{AP,A'^) = A^ - A'P\n[l + A'^/M'^]. From Eq. 
© it is easy to find that when Mq < I/Stt^, this equation 
has two different solutions in the chiral limit. However, 
when m ^ 0, one could only find one solution, the Nambu 
solution, which satisfies M > 0. This result is consistent 
with the one derived from the analysis of the quark DSE 
under rainbow approximation. 

Obviously, the vacuum of the Wigner phase should be 
different from that of the Nambu phase and their differ- 
ence can be characterized by the quark condensate which 
is regarded as the order parameter for chiral phase transi- 
tion. Therefore, the gluon propagator should be different 
due to different values of quark condensate in these two 
phases. In order to refiect this fact, we introduce the 
quark condensate contribution (Eq. ([6|)) to the gluon 
self-energy and modify the effective gluon propagator in 
the normal NJL model as following 
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(10) 

where the first term in the right-hand-side of Eq. (ITU|) is 
the usual model gluon propagator employed in the NJL 
model which has the same form in both Nambu phase 
and Wigner phase and can be regarded as the pure Yang- 
Mills part D^i^^ (q) in the present work; the second term, 
which is inspired by the result of QCD sum rules J12i |. 
is the leading order non-perturbative contribution from 
quark condensate through quark loop insertions. Here, 
it should be noted that according to the usual approx- 
imation of NJL model, in obtaining the current quark 
mass dependent term of Eq. pUj) we have taken all the 



momentum dependence of the effective interaction in mo- 
mentum space as a constant. For this purpose, we have 
introduced a momentum scale A, which reflects the large 
distance behavior of QCD. For the external momentum 
squared much larger than A^ the current quark mass de- 



pendent term of Eq. (|T0|) can be neglected, whereas for 
external momentum squared approaching A^ the contri- 
bution of quark condensate which has been neglected in 
the normal NJL model must be considered. Just as will 
be shown below, the current mass dependent term in Eq. 
(10) plays an important role in searching for the Wigner 
solution at finite current quark mass and the study of 
partial restoration of chiral symmetry at finite chemical 
potential. 

The value of Mq which accounts for the pure Yang- 
Mills gauge field contribution could be fixed by requiring 
the amount of the intensity of the effective interaction to 
be Meff/A = 0.17 for the Nambu solution M = 238 MeV 
which is determined by fitting the observables such as 
pion decay constant and pion mass (in the present paper 
we set the current quark mass m = 5 MeV) [i8[. With 
the modified gluon propagator given by Eq. ([T0|. the 
quark gap equation Eq. 1^ becomes 
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Now let us turn to the the calculation of Eq. (TTTI) . To il- 
lustrate how the solution of Eq. (fTTj) varies with different 



A^, let us define 



F{M) ^ M ~m 
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and the solution of F{M) = is just the solution of the 
quark gap equation (Eq. ([TT|) ). In Fig. [2]F(M) is plotted 
as a function of M with different Kq. From Fig. [2] it can 
be seen that when Kq is larger than about 100 MeV, the 
equation F{M) = has only one solution M = 238 MeV, 
which is similar to the situation discussed in Ref. [3]. 
When A, < 100 MeV, the equation F{M) = has three 
solutions. Specifically, when 70MeV < A, < 100 MeV, 
one solution is the required Nambu solution M = 238 
MeV, and the other two solutions are all smaller than it; 
when Aq < 70MeV, among the two solutions other than 
the Nambu one, one is smaller than it and the other 
one is larger than it. Here we note that physical observ- 
ables require the Nambu solution to be M = 238 MeV 
and the stability condition of the Nambu solution would 
exclude the existence of solutions larger than it. There- 
fore, the parameter A^ should be constrained within the 
range 70MeV < A, < lOOMeV. For A, in this range, the 
smallest solution of Eq. pT|) is not very large compared 
with the current quark mass, and when the current quark 
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FIG. 2. Solutions of the gap equation with different Aq (m is 
fixed to be 5 MeV) 



mass m tends to zero, this solution will continuously tend 
to zero, which is just the Wigner solution in the chiral 
limit. Therefore, this solution might be identified as the 
Wigner solution in the case of m ^ which describes the 
perturbative dressing effect. 

The result in Fig. [5] shows that the scale at which the 
current quark mass dependent term of Eq. (jlOp would 
affect the effective interaction can change the pattern of 
the solutions of the quark gap equation. If the current 
quark mass dependent term plays an important role in 
the infrared region in the effective interaction (Ag < 100 
MeV), then the intensity of the pure Yang- Mills field 
would be weakened and the Wigner solution will appear. 
On the contrary, when A^ > 100 MeV, the pure Yang- 
Mills part would be dominating and strong in the in- 
frared region, and therefore the Wigner solution cannot 
exist due to strong interaction. From physical consider- 
ation we choose A^ = 70 MeV, because in this case the 
gap equation has just two solutions which can be iden- 
tified as the Nambu solution and the Wigner solution. 
For A, = 70 MeV, we plot the F{M) versus M curve 
for different current quark mass m in Fig. [3] It can be 
seen that as m increases, the effective mass M of dressed 
quark in Wigner phase will increase and at last coincide 
with the Nambu solution when m ~ 60 MeV. 



Tabic I The Nambu and Wigner solution 
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As usual, the quark condensate is defined as 
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and its value for the two solutions is listed in Table I (In 



Table I we list the solution of Eq. dTTI) with A^ 
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FIG. 3. Solutions of the gap equation with different m (A, 
fixed to be 70 MeV) 



FIG. 4. Solutions of the gap equation at finite chemical po- 
tential 



MeV and m = 5 MeV) . It can be seen that the value 
of the quark condensate in Nambu phase is larger than 
that of Wigner phase by one order of magnitude, which 
represents DCSB of Nambu phase. Here it should be 
pointed out that the quark condensate of Wigner solution 
is small but non-zero because it reflects the explicit chiral 
symmetry breaking due to non-zero current quark mass. 
Of course, in order to determine which solution is the 
real one, one should compare the pressure (thermody- 
namical potential) of the different solutions. The vacuum 
pressure V of the two solutions are also listed in Table I 
which is calculated via "steepest descent" approximation 
as following Q 



r 



d*p 
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From Table I it can be seen that the vacuum pressure 
of Nambu phase is much larger than that of the Wigner 
phase (more than two orders of magnitude) , which means 
Nambu phase is more stable than the Wigner phase when 
temperature and density are zero. The vacuum pressure 
difference of the two phases can be regarded as the bag 
constant -Bbag and the results in Table I correspond to 
-Bbag ^ (206MeV)'^ which is consistent with the value 
used in the literature Q. One may expect that with 
increasing temperature and/or density this quantity may 
change and chiral phase transition would happen. We 
will discuss this question in the next section. 



IV. PARTIAL RESTORATION OF CHIRAL 

SYMMETRY AT FINITE CHEMICAL 

POTENTIAL 

Now we can generalize the previous treatment to the 
case of finite density. The quark propagator at finite 



quark chemical potential /x could be expressed as follow- 
ing 



G-^{p,^i) = tA^+B-Cfi"f4 



(15) 



where A(pP',p4_, /i), B(jP,p4, /i) and G{p^,p4, /i) are scalar 
functions of p^, p4 and /z. With the model gluon propa- 
gator in Eq. (fTO|) the DSE of quark propagator at finite 
chemical potential is 



lA;^ + B — C/i74 = i;^ + m — 11^4 

4 f d^q 2iA4+AB ~2G^i-i4 



3M; 



ff 



(27r)4 A2g2 -)- ^2 _ ^2^2 ^_ 2iACmi 



,(16) 



From the above equation one could easily find the solu- 
tion should be ^ = 1 and B and C are constant. The 
constant Cfji plays the role of effective chemical poten- 
tial and therefore let us set /z* = C/i and B — M which 
satisfy the following combined equations 
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with quark number density p{lJ,*) defined as following [K 
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p{fi*) = -N,Nf 
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-tr[Giq,fi)^4], (19) 



and EgM = \/F+M2. The Eqs. dH]) and ^ are 
numerically solved and the results are shown in Fig. 
m From Fig. S] One could find the effective mass of 
the dressed quark in the Wigner phase decreases with 




FIG. 5. The Bbag(At) at finite fj. 



increasing fi, which means with increasing density the 
dressing effect of quarks becomes more and more weak. 
On the other hand, the corresponding one in Nambu 
phase decreases with increasing fj, until n ~ 160 MeV, 
and when fi > 160 MeV the effective mass of the dressed 
quark in Nambu phase increases with increasin g ii. 

In the previous literature (see, e.g., Refs. [la,[20|) ^^^ 
usually employ the maximum of the susceptibility gj^' 
to determine the transition temperature. In fact, a more 
reliable criterion for the chiral phase transition is the 
pressure difference of the Nambu phase and the Wigner 
phase, i.e., the bag constant Bha,g{p). The pressure den- 
sity of the two solution at finite chemical potential could 
be calculated as following [Tg 



v{fi) = p(^ = 0) + / d^l'p{fl') 



(20) 



where the pressure density of the vacuum 7'(/i = 0) can 
be calculated through Eq. (ITU) . The i3bag(A*) is plotted in 
Fig. [S]in which one can see when fi < fj,c = 260 MeV the 
Nambu solution is more stable and when ^ > fic = 260 
MeV the Wigner solution is more stable. At Hc — 260 
MeV the pressure of the two phases is equal and the two 
phases could co-exist at this point. Here it should be 
noted that no one has calculated the i?bag(M) in the case 
of nonzero current quark mass in the past. This is due 
to lack of knowledge about the Wigner solution of the 
quark gap equation at finite current quark mass in the 
previous literature. 

Here it is interesting to compare our results with those 
of the normal NJL model. The first-order phase tran- 
sition point /Zc in our modified NJL model is smaller 
than the one obtained in the normal NJL model which 
is about 354 MeV or 500 MeV corresponding to different 
parameters [2l| . It should also be pointed out that in the 
normal NJL model, the second solution appears when fi 
is big enough 18], but the magnitude of this solution 
at the phase transition point is much bigger (about 110 
MeV or 130 MeV, see Ref. ^]) than the Wigner solu- 



tion obtained in the present paper (about 15 MeV). In 
addition, we want to stress that the Wigner solution at 
finite current quark mass in the normal NJL model is 
due to density effect. When the chemical potential tends 
to zero, this solution disappears. This shows fi — is a 
singularity of the Wigner solution at finite current quark 
mass. If this is real, it means that one cannot study the 
Wigner solution by means of small /i expansion, while 
the method of small /x expansion is a usually employed 
one in the study of lattice QCD at finite density. 



V. SUMMARY 

To summarize, based on the general analysis that the 
dressed gluon propagator in Wigner phase should be dif- 
ferent from that in Nambu phase, we introduce the con- 
tribution of quark condensate to the gluon propagator 
and investigate the solution of quark DSE in the case of 
nonzero current quark mass. With such a modified model 
gluon propagator, in the framework of NJL model we 
show that the quark DSE indeed has a Wigner solution in 
the case of nonzero current quark mass. We then gener- 
alize this approach to the case of finite chemical potential 
and discuss partial restoration of chiral symmetry at fi- 
nite chemical potential. From the calculated result of the 
bag constant we find that when /x < /Ltc = 260 MeV the 
Nambu solution is more stable and when jjl > jjLc — 260 
MeV the Wigner solution is more stable. At fic — 260 
MeV the pressure of the two phases is equal and the 
two phases could co-exist at this point. We also com- 
pare our results with those of the normal NJL model. It 
is found that the first-order phase transition point ^c in 
our modified NJL model is smaller than the one obtained 
in the normal NJL model which is about 354 MeV or 500 
MeV corresponding to different parameters. In addition, 
in the normal NJL model, the second solution appears 
when /i is big enough, but the magnitude of this solution 
at the phase transition point is much bigger (about 110 
MeV or 130 MeV) than the Wigner solution obtained in 
the present paper (about 15 MeV). Finally we would like 
to point out that the results obtained in this paper are 
based on a simple NJL model. It is well-known that the 
NJL model is far from QCD. In order to obtain a more 
solid result, we need to further discuss this problem in 
the framework of a model with better QCD foundation, 
such as DSE of QCD 0. 
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